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Abstract 

The variational principle for linear stability of three-dimensional, in- 
homogenious, compressible, moving magnetized plasma is suggested. The 
principle is "softer" (easier to be satisfied) than all previously known varia- 
tional stability conditions. The key point of the analysis is a conservation 
in variations of new integrals inherent in the linearized equation of the 
motion that was not earlier discussed in the literature. 

PACS46.15.Cc, 52.30.Cv 

It is well known that stability of the static equilibrium of magnetized plasma 
can be described by so-called "energy principle" [1]. The principle claims that if 
the second variation of the potential energy, W, of the system "plasma-magnetic 
field" is positive definite near the equilibrium point, then this point is stable. 
Sufficiency of the claim follows from the Lyapunov stability theorem, and neces- 
sity can also be proved [1] . Note that the above mentioned second variation of 
the potential energy corresponds exactly to the potential energy of the linearized 
equation system. 

The main drawback of the principle is that there always may be neutral 
perturbations which do not perturb any physical quantity - and, therefore, W 
as well. Thus, the second variation of the potential energy can be guaranteed 
to be only positive semi-definite. In other words, using the energy principle [1], 
one can talk about spectral stability only, namely, about presence or absence of 
imaginary frequencies in the spectrum of the linearized force operator (nonlinear 
stability needs an analysis of neutral perturbations - see, e.g., [2]). Continuing 
this logic, we restrict ourselves with linearized equations. 

The attempt of using the similar approach to investigate stability of moving 
plasma performed by Frieman and Rotenberg [3] was not so lucky, although the 
energy principle was formally obtained. Their result can be briefly described as 
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follows. Consider the linearized equation of motion for plasma displacement £ 
in the frame of ideal one-fluid magnetohydrodynamics, 



where the linearized force operator, 

F(£) = -<Jp(V-V)V - i o(<5V-V)V - p(V-V)5V 
- VSp + (V x SB) x B + (V x B) x SB , 

is combined of usual perturbed quantities, 

s P = -vipti), sv = (v-v)€ - (*-v)v , 

5p = Vp - ^p V-£ , (5B = V x x B) 

(note that 5V denotes here only the part of full Eulerian velocity perturbation - 
the part, which survives even for time- independent displacements, £). Station- 
ary plasma density, p, velocity, V, pressure, p, and magnetic field, B, satisfy 
the following equilibrium conditions: 



Dot means a partial time-derivative, 7 means the adiabatic exponent. Force 
operator is proved to be self-adjoint in the following sense, 



p£ + 2p(V-V)£-F(£) = 



(1) 



p(V-V)V + Vp = (V x B) x B , 

V-(pV) = , 

V- Vp + 7p V V = , 

V x (V x B) = . 




while the second term in (1) is obviously antisymmetric: 




Multiplying Eq. (1) by £ and integrating over the whole space, we found the 
energy conservation in the form E — 0, where 




Minimizing E over we approach to the energy principle by Frieman-Rotenberg, 




(3) 
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Contrary to the static case (V = 0), in which condition (3) appears to be both 
sufficient and necessary for linear stability in the case of V ^ 0, condition (3) 
is normally too strong, and never can be satisfied except of field-aligned flows 
(V ~ B)[3] or of those which may be reduced to the field- aligned flows (see, 
e.g., W). 

Energy principle (3) may be improved by use of the Arnold conjecture [5]- 
[7], following which we have to add to the energy (2) the set of other known 
integrals of the motion. Speaking in other words, variables £ and £ in (2) 
are not absolutely independent but subject to the constraints resulting from 
conservation of other integrals of the motion. 

Such an improved principle was derived by Ilgisonis and Pastukhov [8] , then 
it was verified by Hameiri [9]. It was also re-obtained [10] with help of Pfirsch & 
Morrison's method [11] of dynamically accessed perturbations. That stability 
condition is currently the best among the known ones, although it is still not 
appropriate for arbitrary stationary plasma flow. 

For the linearized equation (1), the Ilgisonis & Pastukhov extra invariants 
can be written in terms of neutral perturbation £ N : 

F(€ JV ) = 0, d t £ N = o. 

Multiplying Eq. (1) by £ N and integrating again over the space, we have 1 = 0, 
where 

I = J {p'Hn + 2p£iv(V-V)£) d\ . (4) 

For the system with nested set of magnetic surfaces, -0(r) = const, £ N may be 
generally represented as 

€n = A u (V0u + A„(V0v , (5) 

where u = B/p, v = D/p, and D is a divergence-free frozen- in-plasma vector, 
tangential to the same magnetic surfaces, -0(r) = const, but different from B, 

B x D = pVip (6) 

- see [8, 9] for explanations of how D can be built-up. For = in (5), con- 
servation of / corresponds to the cross-helicity invariance, which is well known 
contrary to more general quantity /. Note that taking into account (5), (6), the 
second term under the integral in Eq. (4) can also be written as 

or as 

/^•(V-V^ + pV-tov-V)*, 

for absolutely arbitrary functions A Ujt , (tp). 
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It is very important that Eq. (1) allows for some extra set of invariants 
different from (2) and (4). Indeed, differentiating Eq. (1) with respect to time, 
then multiplying it by £ and integrating, we found - like in the case of energy but 
using once again original equation of motion (1) - that the following quantity, 

E i = \j {^(£)-2pV-V£) 2 -£.F(£)j d\ , (7) 

is conserved. This invariant is exact for linearized dynamics (1), and cannot 
be reduced to the conservation of energy (2). In principle, we may continue 
the procedure and get in the same manner an infinite set of similar invariants. 
However, to investigate a stability, it might be sufficient to involve into our 
analysis only finite number of the invariants. Here we show that taking into 
account even one of them, E 2 , it appears to be possible to improve the stability 
condition significantly. 

Note that being varied separately (i.e., when £ is independent), invariant (7) 
results in the same Fricman-Rotenberg condition of semi-positive definiteness 
of quadratic form (3) based on F-operator. To improve stability condition, we 
will consider the functional 

U = E + pE 2 - I (8) 

to be varied over £ and £, subject to the independent conservation of the inte- 
grals, E 2 and /. Explicitly, 

°-f {f ? + |h-^)'-f«, 

+(2p£-(V-V) - P i)(Ku + A„v) j d 3 r , (9) 

where the constant, ^, and 1-D functions, X UiV (il)), play roles of the Lagrangian 
multipliers; we have to choose them to provide the integrals E 2 , I be equal to 
their equilibrium values, i.e., to zero. 
Functional U is minimized by £: 

£ = A„u + A„v + 2pV-V ( 2(V-V)£ - . (10) 

v . ' p V p J 

is 

Putting /i — > in Eqs. (9), (10), we approach to the Ilgisonis-Pastukhov-Hameiri 
condition [8, 9]. Indeed, we have at fi — > 0: 
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where £ N satisfies (5), and \ u ,v'- 

(^-2pe(V-VK N )rf 3 r = 0. 



/< 



In this limit, we have not really used the condition of ^-conservation. Note that 
Uiph > — J d 3 r/2, and, therefore, the Ilgisonis-Pastukhov-Hameiri con- 

dition, Uiph > 0, is "softer" than the condition (3) by Frieman and Rotenberg. 
However, it is still not appropriate for arbitrary flow, in which V is not parallel 
to B. The sign-indefinite term, £-F(£) , contains the high-order ^-derivatives 
and, therefore, can always prevail on the positive term p(? N (see [9] for more 
details). 

Now let us account for small but finite /x. Solving Eq. (10) by iterations in 
ti, we found 

C » €jv - 2/i(V-V)€ , (11) 

where we used the notation 

^0 = F(€)/P-2(V-V)€ JV . (12) 
Stability condition is expressed again by functional U depending on £: 

U = j {f {t N 2 At (V.V)^ ) 2 - d 3 r>0. (13) 

Here £ N and £ are defined by Eqs. (5), (12), and also depend on 

Lagrangian multipliers have to be found by substituting Eq. (11) into the 
conditions 



o 



4=€ w -2m(v-V)£ 



They are: 



r 



_ A,£> M - A D V _ A U D V - A £> M 

u ~~ 4 4 - A 2 ' v ~ A A - A 2 ■ [ ' 



Here 



4=u,. =< 4itp((V-V)w) 2 - pw 2 > , 

A Q =< 4/xp((V-V)u)-((V-V)v) - pu-v > , 

=< 2(a*F(€) - p£HV-V)w > , 

and angular brackets mean the averaging over magnetic surface. 

Note that the left-hand-side of the condition (13) contains the high-order 
derivatives of £ in the first (non-negative) term, hence, the second (sign-indefinite) 
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term is no more critical. It is the main advantage of the condition (13) with 
respect to previous one, Uiph > 0, that it has a sense for arbitrary (not only 
field-aligned) flow and, therefore, may have a practical merit. 

This work was partially supported by the Human Capital Foundation Grant 
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